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Abstract: 

All lowest-order amplitudes for e"''e~ — > 4/7 are calculated including five anomalous 
quartic gauge-boson couplings that are allowed by electromagnetic gauge invariance and 
the custodial SU(2)c symmetry. Three of these anomalous couplings correspond to the 
operators £o; and £„ that have been constrained by the LEP collaborations in WW7 
production. The anomalous couplings are incorporated in the Monte Carlo generator 
Racoon WW.0 Moreover, for the processes e"'"e~ — 4/7 Racoon WW is improved upon 
including leading universal electroweak corrections such as initial-state radiation. The 
discussion of numerical results illustrates the size of the leading corrections as well as the 
impact of the anomalous quartic couplings for LEP2 energies and at 500 GeV. 
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1 Introduction 



In recent years, the experiments at LEP and tlie Tevatron liave establislied tlie ex- 
istence of elementary self-interactions among three electroweak gauge bosons, mainly by 
analysing the reactions e+e~ W+W^ and pp — » W7 + X. The empirical bounds (see 
e.g. Ref. [ |1|) on anomalous triple gauge-boson couplings confirm the Standard- Mo del 
(SM) couplings at the level of a few per cent. Recently, the LEP collaborations have 
started to put also bounds on anomalous quartic gauge-boson couplings (AQGC) upon 
studying the processes e"'"e~ — >■ W"'"W~7, e"'"e~ — > Z77, and e"^e~ —>■ uu'-f'-f. The OPAL 
[ H , L3 [ ^ and ALEPH collaborations have already presented first results, which have 
been combined by the LEPEWWG [ |. 

The experimental analysis of anomalous triple and quartic gauge-boson couplings re- 
quires precise predictions from Monte Carlo generators including these anomalous cou- 
plings. In particular, it is necessary to account for the instability of the produced weak 
bosons, which decay into fermion-antifermion pairs. While several generators including 
triple gauge-boson couplings in e^e~ WW — > 4/ exist for quite a long time [ || , up to 
now no generator has been available that deals with the processes e+e~ —>■ WW(7) —>■ 4/7 
in the presence of AQGC.Q As a preliminary solution [ |^, Q a reweighting technique was 
used in existing programs for e"'"e~ —>■ WW(7) —>■ 4/7 in the SM where the SM matrix 
elements were reweighted with the anomalous effects deduced from the program EEWWG 
[ 1^ for on-shell WW7 production. The aim of this paper is to improve on this situation 
by implementing the relevant AQGC in the Monte Carlo generator Racoon WW [ ^, 
which is at present the only generator for all processes e"'"e~ — > 4/7. SM predictions 
for all 4/7 final states obtained with this generator were presented in Ref. [ 0; further 
results for specific final states can be found in Refs. [ H, |10|, |ll|. Here we supplement 
these numerical results by a study of AQGC effects at LEP2 and linear collider energies. 
Moreover, we include two quartic gauge-boson operators in the analysis that, to the best 
of our knowledge, have not yet been considered in the literature before. 

As a second topic, we improve the RacoonWW predictions for 4/7 production by in- 
cluding the dominant leading electroweak corrections. In particular, we take into account 
additional initial-state radiation (ISR) at the leading-logarithmic level in the structure- 
function approach of Ref. [ , where soft-photon effects are exponentiated and collinear 
logarithms are included up to order 0{a^). Leading universal effects originating from the 
renormalization of the electroweak couplings are included by using the so-called G^-input- 
parameter scheme. The singular part [ |T3| of the Coulomb correction, which is relevant for 
intermediate W-boson pairs in e"'"e~ — > WW(7) 4/7 near their kinematical threshold, 
is also taken into account. 

The paper is organized as follows. In Section ^ we introduce the relevant AQGC 
and give the corresponding Feynman rules, which are used in Section ^ to calculate the 
AQGC contributions to all e"'"e~ — > 4/7 amplitudes. In Section ^ we improve the tree-level 
predictions by including leading universal electroweak corrections. Section ^ contains our 
discussion of numerical results, which illustrate the impact of the leading corrections to 
the SM predictions as well as the effects of the AQGC. A summary is given in Section 1^. 



-'^ While finishing this paper a version of the Monte Carlo generator WRAP [ g] became available that 
also includes AQGC but uses a different set of operators. 
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2 Anomalous quartic gauge-boson couplings 



Since we consider the class of e+e~ 4/7 processes in this paper, we restrict our anal- 
ysis to anomalous quartic gauge-boson couplings (AQGC) that involve at least one photon. 
Moreover, we consider only genuine AQGC, i.e. we omit all operators that contribute also 
to triple gauge-boson couplings, such as the quadrilinear part of the well-known opera- 



tor F'^'^W^'^Wp^. Imposing in addition a custodial SU(2)c invariance [ to keep the 
p parameter close to 1, we are left with operators of dimension 6 or higher. Following 
Refs. [ 1^, [16|, [l^ we consider dimension-6 operators for genuine AQGC that respect 



local U(l)cm invariance and global custodial SU(2)c invariance. These symmetries reduce 
the set of such operators to a phenomelogically accessible basis. More general AQGC 
were discussed in Ref. \ \l 



In order to construct the relevant AQGC, it is convenient to introduce the triplet of 
massive gauge bosons 

W, = {Wl,Wl,Wl) = (^-^{W^ + W~),, -^{W - W-),, ^Z,^ , (2.1) 

where W^j^ and are the fields of the and Z bosons, and the (abelian) field-strength 
tensors 



W'"'' = d^'W''' - d'^W''", (2.2) 

where A^ is the photon field. The parameter c„ is the cosine of the electroweak mixing 
angle. The quartic dimension-6 operators are obtained upon contracting two factors of 
with two field-strength tensors. Under the explained symmetry assumptions there 
are five independent AQGC operators of dimension 6. We choose the following basis: 

Co = -^aoF'^'^F^^WjW", 
Cc = -^a^F^^F^pW^W^, 

4 = - e,,, F^^WlJViW'''', (2.3) 



F^, = -e^.p^F"^ (£°i2^ = +1) (2.4) 



where 

1 

— L 

2 

is the dual electromagnetic field-strength tensor, and e is the electromagnetic coupling. 
The scale A is introduced to keep the coupling constants Oj dimensionless. The operators 
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£o and which were introduced in Ref. [15], conserve the discrete symmetries^ C, P, 
and CP, while the others respect only one of these symmetries. The operator which 
was defined in Refs. [ 0, |16|, |T^, conserves only P, but violates C and CP. The P- violating 



operators Cq and £n have to our knowledge not yet been considered in the literature. 
While Cq conserves C and violates CP, £n conserves CP and violates C. 

We add some comments on the completeness of the set ( |2.3|) of quartic couplings. 
At first sight, there are three more P- violating couplings of dimension 6 that can be 
constructed with the tensor e^^p^, namely 

(2-5) 

These operators can be reduced to Cq and £n by exploiting the Schouten identity 

(2.6) 

which is a consequence of the four-dimensionality of space-time. Moreover, we could 
have constructed also operators from d^Wi, and d^^V^ separately instead of taking W^,^. 
However, the new operators obtained this way only lead to additional terms involving 
either d^W^ or d^F^y, which do not contribute to the amplitudes for e"'"e~ — > 4/7 for 
massless external fermions. 

In order to deduce the Feynman rules for the considered AQGC, we express the fields 
W^^ in terms of physical fields. 



K^^iw-Zp - w-z^) - w-yiw^Zp - w;z^) 

+ Z,yiW+W~-W;W-)]. (2.7) 



Taking all fields and momenta as incoming, we obtain the Feynman rules 
{W+;Z^},p3 





x{Aao9''''[iPiP2)9'"'-PiP'2] 

+ ac Mp^2 + pfP2)9"' + {P1P2) {9''' 9"^ + 9"'' 9'^) 



^We adopt the usual convention [ |9| that PVfj,P^^ = and CV^C^^ = —V^ for aU electroweak 
gauge bosons. While for the photon {V — A) these transformations follow from the C and P invariance of 
the electromagnetic interaction, for the weak bosons {V — W^,Z) they are mere definitions, which are, 
however, in agreement with the C and P invariance of the bosonic part of the electroweak interaction. 
The CP transformation, on the other hand, is well-defined for all electroweak gauge bosons. 
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+ ^~aog''^PipP2aS'''n, (2-8) 




16A2c, 



Zu,P2 Wp,p4 

x{a„ [-(piP2)(^?'^°/^ - g'^gn - {piPsKg'^g"'' - g^'g'^^) - {prnKg^'g''^ - g'^^g^n 
+ pWg^" - pfgn + pMg'''' - plg""^) + p'Mg""^ - Pig^l 

- g'^iP^iP^s - P'ipl) - g^'-ip'ipl - pfe) - g'^W2 - p'ipt)] 

+ ~anPiA{Pi+P2re''^^' + {pi+PzTe^'^P + (^1+^4)^"'^'' 

- {P2 - PsW'^s^^^' - iPs - PaW^s^'-^p - (p4 - P2)./"£™^1}. (2.9) 

Note that the 'yZW^W" couphng is symmetric with respect to cychc permutations of 
ZW+W-/i.e. of (p2,^^),(P3,«),(P4,/3). 

In order to evaluate the diagrams with the P-violating couphngs within the Weyl-van 



der Waerden spinor formahsm (see Ref. [ |20| and references therein), which we use in 



the calculation of our amplitudes, the tensor e'^'^p'^ has to be translated into the spinor 
technique. Following the notation of Ref. [ ^ the tensor is substituted in the Feynman 
rules according to the identity (e°^^^ = +1) 

e'^'^P'^ [la^^) [la^^) [la^^) [la^^) = 1^(^^ae ^cg ^bd ^fh _ ^ac ^eg ^bf ^Duy (2.IO) 

For the the Standard-Model (SM) parameters and fields, i.e. for the SM Feynman 
rules, we follow the conventions of Ref. [ ^ 

3 Amplitudes with anomalous quartic gauge-boson couplings 

In Ref. [ ^ we have presented the SM amplitudes for all e"'"e~ 4/7 processes with 
massless fermions in a generic way. The various channels have been classified into charged- 
current (CC), neutral-current (NC), and mixed CC/NC reactions, and all amplitudes 
have been generated from the matrix elements A^cca and TVlNCa, which correspond to 
the simplest CC and NC reactions, called CCa and NCa, respectively: 

CCa: e+e- ^ //'FF', 

NCa: e+e- //FF, 



•^In this context it is important to recall that different conventions are used in the literature concerning 
the sign of the electroweak gauge coupling g = e/sw, and thus of the sign of the sine of the weak mixing 
angle Sw Since the SM quartic coupling 7ZW+W- changes under the inversion of this sign, which of 
course can never affect physical quantities, the anomalous jZW^W^ coupling also has to be reversed 



when switching from one convention to the other, although does not appear in (2.9) explicitly. 
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fa{Pa, (Ja) 




fb{Pb, CTb 



Figure 1: Generic diagram for the AQGC contribution to e"'"e — > 4/7 

where / and F are different fermions (/ 7^ F) that are neither electrons nor electron 
neutrinos {f,F 7^ e~,z/e) and their weak-isospin partners are denoted by /' and F', re- 
spectively. Here we supplement the SM amplitudes of Ref. [ ^ by the corresponding 
contributions resulting from the AQGC given in ( p.3|) . We follow entirely the conventions 
of Ref. [ and denote the external particles of the considered reaction according to 

e+(p+, cr+) + e"(p_, o-_) fi{ki, ai) + f2{k2, as) + Mh, as) + Uih, 0-4) + 7(/c5, A), 

(3.1) 

where the momenta and helicities are given in parentheses. We list the expressions for the 
contributions 5A^cca,AQGC and 5A^NCa,AQGC of the anomalous couplings to the generic 
CC and NC matrix elements A^cca and A^NCa, respectively, which have to be added to 
the SM contributions. From A^cca and AlNCa the amplitudes for all other CC, NC, and 
CC/NC reactions are constructed as explained in Ref. [|^. Some minor corrections to 
this generic construction are given in App. 0. 

We express the AQGC contributions 5A^cca,AQGC and 5AiNCa,AQGC in terms of the 
two generic functions Al7yy,AQGC and A4zww,acigc, which correspond to the •yyVV and 
■jZWW couplings, respectively, with V = W, Z, 

'5-^ccZaqgc '"''"''^(P+'P-' ^2, h, h) 

= A^^wXqgc'"''"'''' ~^2, -h, -h, h) 

+ -^zw,AQGc''''"''''"''*'^(P+'P-' -^1' -^2, -h, -h, h), (3.2) 

^-^NCa,AQGC '""'""'^(P+'P-' h, ^^2, h, h, h) 

= ^7ZZ,AQGC~'''' '^(P+'P-' ~^2, -fcs, -^4, h) 

+ '^7ZZ,AQGc'''~''''''^'''"'^(-^l> -^2, -h, -k4,P+,P-, k^) 

+ A^^IzJqgc "^"'"''""'^"'(-^3, -h,p+,p-, -fci, -k2, k,). (3.3) 

The generic Feynman graph that corresponds to A^yiVaVa.AQGC is shown in Figure where 
the fermions and antifermions are assumed as incoming and the photon as outgoing. 
Explicitly the generic functions read 

•^7yy!AQGc''"'''^'^(Pa, Pb, Pc Pd, Pe, Pf, k) 

^ ~ ^jf^ ^^^^~-^^-^^-'^^-^'--f 9^lhf,9v%fyVf.ffPviPc + Pd)Pv{Pe + Pf) 



X 



8A 



{Pa,Pb,Pc,Pd,Pe,Pf,k) + 



{Pa,Pb,Pc,Pd,Pe,Pf,k) 



with Cj^fzz = 



■yyWW 



(3.4) 



M 



0'a,0-fj,(Tc,(Jd,Ue,0-f,\ 

ZWW,AQGC 



iPa,Pb,Pc,Pd,Pe,Pf,k) 



le 



X PziPa + Pb)Pw{Pc + Pd)Pw{Pe + Pf) A''^l'''^'''^'\pa,Pb,Pc,Pd,Pe,Pf,k), (3.5) 

where the propagator functions Pvip) and the fermion-gauge-boson couphngs gyjfi can 
be found in Ref. [ |^. We have evaluated the auxihary functions ^4^^'°"'=''^''''^ with = 0, c, n 
in terms of Weyl-van der Waerden spinor products [ EH]: 



At^^^{Pa,Pb,Pc,Pd,Pe,Pf,k) = (ao + iflo) 



{{Pbky f {PdPf)*{PcPe) 
(PaPb)* 



(3.6) 



'{Pa,Pb,Pc,Pd,Pe,Pf,k) 



(Pa ■ Pb) 

X \^{PaPbY{Pdk)*{pfkY{paPc){PaPe) + {{PbkYf {PdPf)* {PaPb) {PcPe)], (3.7) 

^t^'^'^iPa,Pb,Pc,Pd,Pe,Pf,k) = (ttn + iSn) 

X { {PbPf)*{Pdk)*{PaPe) {Pak)*{PaPc) + {Pbk)* {pbPc) + {Pek)* {PcPe) + {Pfk)*{PcPf) 

+ {PbPd)*{Pfk)*{PaPc) {Pck)*{PcPe) + {Pdk)*{pdPe) " (Pak)* (PaPe) " (Pbk)* (pbPe)] 

- {PdPf)*{Pbk)*{PcPe)[{Pek)*{PaPe) + {Pfk)*{PaPf) - {Pck)* {PaPc) - {Pdk)* {PaPd)]} ■ 

(3.8) 

The remaining polarization combinations follow from crossing and discrete symmetries, 





^{Pa,Pb,Pc, 


Pd,Pe, 


Pf^ 


k) 


= 

o-k 


'{Pa,Pb,Pc,Pd,Pf,Pe,k), 


^~^iPa,Pb,Pc, 


Pd,Pe, 


Pf^ 


k) 


= At^^' 

o-k 


'iPa,Pb,Pd,Pc,Pe,Pf,k), 




^{Pa,Pb,Pc, 


Pd,Pe, 


Pf, 


k) 


= K^^' 

a-k 


'{Pa,Pb,Pd,Pc,Pf,Pe,k), 




^{Pa,Pb,Pc, 


Pd,Pe, 


Pf^ 


k) 




'{Pb,Pa,Pc,Pd,Pe,Pf,k), 


O-k 


^iPa,Pb,Pc, 


Pd,Pe, 


Pf^ 


k) 


o-k 


'iPb,Pa,Pc,Pd,Pf,Pe,k), 


\k ^^iPa,Pb,Pc, 


Pd,Pe, 


Pf, 


k) 


= K^^' 

a-k 


'{Pb,Pa,Pd,Pc,Pe,Pf,k), 




^{Pa,Pb,Pc, 


Pd,Pe, 


Pf^ 


k) 




'{Pb,Pa,Pd,Pc,Pf,Pe,k), 




~iPa,Pb,Pc, 


Pd,Pe, 


Pf^ 


k) 




~''^'~'''''^iPa,Pb,Pc,Pd,P 



, A; = 0, c, n. 

(3.9) 

It is interesting to observe that the helicity amplitudes for oq and oq, and similarly for 
and On, differ only in factors ±i for equal coupling factors. These AQGC are the ones that 
are related by interchanging a field-strength tensor F with a dual field-strength tensor F 
in the corresponding operators in ( p.3|) . 

As we had already done in Ref. [ |^ in the case of the SM amplitudes, we have numer- 
ically checked the amplitudes with the AQGC against an evaluation by Madgraph [ , 
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which we have extended by the anomalous couphngs. We find numerical agreement for a 
set of representative 4/7 final states. 



4 Leading universal electroweak corrections 

Besides the genuine AQGC we have also included the dominant leading electroweak 
corrections to e"'"e~ —>■ 4/7 into RacoonWW, similar to our construction [ ^ of an 
improved Born approximation (IB A) for e"'"e~ WW — > 4/. 

The dominant universal effects originating from the renormalization of the electroweak 
couplings are included by using the so-called G^-input-parameter scheme. To this end, 
the global factor in the cross section is replaced by aQ^a{0) with 

o:a = (4.1) 

While the fine-structure constant a(0) yields the correct coupling for the external on-shell 
photon, ttG^ takes into account the running of the electromagnetic coupling from zero to 
Myif and the leading universal mt-dependent corrections to CC processes correctly. The 
mt-dependent correction to NC processes are not included completely. These could be 
accounted for by introducing an appropriate effective weak mixing angle. However, we 
prefer to keep the weak mixing angle fixed by = 1 — = M^/M|, in order to avoid 
potential problems with gauge invariance which may result by violating this condition. 

Initial-state radiation (ISR) to e"'"e~ — > 4/7 is implemented at the leading-logarithmic 
level in the structure-function approach of Ref. [ |12| as described for e"'"e~ —>■ 4/ in 
Ref. [ I in equations (5.1)-(5.4), 

lda^Bl~^'^' = £dx, |^'dx2r^L(xi,Q2)rLL(x2,Q2) J da^^l^'^'' {x,p+, x,p.). (4.2) 

In the structure function T]:^{x,Q'^) [ H, ^ soft-photon effects are exponentiated and 
collinear logarithms are included up to order 0{a^). The QED splitting scale is a free 
parameter in leading-logarithmic approximation and has to be set to a typical momentum 
scale of the process. It is fixed as = s by default but can be changed to any other scale 
in order to adjust the IBA to the full correction or to estimate the intrinsic uncertainty 
of the IBA by choosing different values for Q^. 

For processes with intermediate W-boson pairs, e'^e" — >■ WW(7) 4/7, the singular 
part [ of the Coulomb correction is taken into account, i.e. in this case we have 

dcT^BA^'^' = da^orn'^'^^ [1 + ^Coui(/, e_)g{P)] , s' = (A;+ + k^)'. (4.3) 

The Coulomb singularity arises from diagrams where a soft photon is exchanged between 
two nearly on-shell W bosons close to their kinematical production threshold and results 



in a simple factor that depends on the momenta k± of the W bosons [ [T^, ^ 



Ocn^AS , ftl, fcl = -yr- Im < In = — 



/9 



s'2 + ki + kt- 2s'kl - 2s'kl - 2klk'^ 
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Figure 2: Generic diagrams contributing to the Coulomb singularity in e^e 4/7 



state (see Figure ^a) or from the final state (see Figure ^b) contribute to the Coulomb 
singularity. Therefore, it is not just given by a factor to the complete matrix element. 
However, applying different correction factors to different diagrams would violate gauge 
invariance. Therefore, we decided to use an effective treatment that takes into account 
the dominant effects of the Coulomb singularity. We actually implemented two different 
variants: 

1 . In the first variant we multiply the complete matrix element with the Coulomb cor- 
rection factor with = ki + k2 and /c_ = k^ + k^. In this way we multiply the 
correct Coulomb correction to all diagrams with ISR (Figure]^). However, in this 
approach we do not treat the Coulomb singularity in diagrams with final-state radi- 
ation (Figure 0d) properly. Nevertheless, this recipe should yield a good description 
of the Coulomb singularity, since the diagrams with two resonant W bosons and 
photon emission from the initial state dominate the cross section. This expectation 
is confirmed by the numerical results presented below. 

2. In the second variant we improve on this prescription by differentiating between 
initial-state and final-state radiation according to the invariant masses in the final 
state. To this end, the W-boson momenta entering the Coulomb correction factor 
are fixed as 




A 



s' 



(4.4) 



This correction factor is multiplied with the auxiliary function 




(fci + k2,k3 + ki) 

{ki + k2 + /cs, k^ + k4) 



for A12 < A125, A34 < A345, 
for A12 > A125, A34 < A345 or 



(fc+,A;_) 



A12 > A125, A34 > A345, A125 < A345, 
for A12 < A125, A34 > A345 or 



{ki + k2, k^ + kA + /cs) 



A12 > A125, A34 > A345, A125 > A345, 



(4.6) 
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where Aj^ = \{ki + kj)'^ — M^\ and Aj^; = \{ki + kj + ki^ — M^\. In this way we 
effectively apply the correct Coulomb correction factor to all dominating doubly- 
resonant contributions, shown in Figure 0. 

The Coulomb singularity is not included in processes that do not involve diagrams with 
two resonant W bosons. 

Finally, we optionally include the naive QCD correction factors (1 + as/ir) for each 
hadronically decaying W boson. 

In order to avoid any kind of mismatch with the decay, Fw is calculated in lowest 
order using the scheme. This choice guarantees that the "effective branching ratios" , 
which result after integrating out the decay parts, add up to one when summing over 
all channels. Of course, if naive QCD corrections are taken into account, these are also 
included in the calculation of the total W-boson width. 



5 Numerical results 



For our numerical analysis we take the same SM input parameters as in Refs. [ P, [TT[ . 
We use the constant-width scheme, which has been shown to be practically equivalent to 
the complex- mass scheme for the considered processes in Ref. [ |^. The errorbars shown in 
the plots for the relative corrections result from the statistical errors of the Monte-Carlo 
integration. 



5.1 Comparison with existing results 

We first compare our results for e"''e~ —>■ 4/7 including leading corrections with results 
existing in the literature. 

Predictions for e"'"e~ 4/7 including ISR corrections have been provided with the 
program WRAP [ 0. First results have been published in Ref. [ |1T| where also a com- 



parison with RacoonWW at tree level was performed. Here we present a comparison 
between WRAP and RacoonWW for the same set of input parameters and cuts as in 



Section 5.2. of Ref. [ ^ but including ISR. In this tuned comparison, the W-boson width 
is kept fixed at Fw = 2.04277 GeV, and neither the Coulomb singularity nor naive QCD 
corrections are included. The results are shown in Figures ^ and ^. The absolute pre- 
dictions on the left-hand sides are hardly distinguishable. The relative deviations shown 
on the right-hand sides reveal that the agreement between WRAP and RacoonWW 
is at the level of the statistical error of about 0.2%. The comparison has been made for 
collinear structure functions. Unlike px-dependent structure functions, coUinear struc- 
ture functions do not allow to take into account the Bose symmetry of the final-state 



photons resulting in some double counting [ However, for not too small cuts on the 



photon energy and angle these effects are beyond the accuracy of the leading-logarithmic 
approximation. This has been confirmed by the numerical analysis in Ref. [ P] . 

In Figures |^-|^ we repeat the comparison between YFSWW3-1.14 (scheme A) [ ^ 



and RacoonWW given in Section 4.1. of Ref. [ |TT| for the photonic distributions. But 
now we include besides the tree-level predictions of RacoonWW for e+e" 4/7 also 
those including leading-logarithmic ISR. Note that unlike in all other distributions dis- 
cussed here, a recombination of photons with fermions is performed for this comparison. 
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Figure 3: Total cross section including leading-logarithmic ISR corrections for the process 
e"'"e~ ud/i~z/^7 as a function of the CM energy with minimal energy of the observed 
photon of 1 GeV. Absolute predictions from WRAP [ ^ and RacoonWW are shown 
on the left-hand side, the relative differences between the two programs are shown on the 
right-hand side. 



a[fb] 
70 



WRAP 



65 - 
60 
55 - 
50 - 
45 - 
40 - 
35 
30 



1 



1 1 

RacoonWW I — \ — I 

WRAP '-^-y 



_i_ 



[GeV] 




[GeV] 



Figure 4: Total cross section including leading-logarithmic ISR corrections for the process 
e"'"e~ — >■ ud/i~z>^7 as a function of the minimal energy of the observed photon for y/s = 
192 GeV. Absolute predictions from WRAP [ H and RacoonWW are shown on the left- 
hand side, the relative differences between the two programs are shown on the right-hand 
side. 
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Figure 5: Distribution in the photon energy for the process e^e^ —>■ ud/i^z/^7 at y/s = 
200 GeV 



We restrict ourselves to the "bare" recombination scheme (see Refs. [ |TT| for details). 
Moreover, the W-boson width is calculated including the full 0{a) electroweak corrections 
together with naive QCD corrections resulting in Fw = 2.08699 GeV. We compare the 
distributions in the photon energy E^, in the cosine of the polar angle 6^ of the photon 
w.r.t. the e"*" beam, and in the angle 6yf between the photon and the nearest charged 
final-state fermion for the process e"'"e~ — > ud/i~z/^7 at y/s = 200 GeV. The differences of 
15-20% between YFSWW3 and the pure Born prediction of RacoonWW (RacoonWW 
Born), which have already been shown in Ref. [|ll|, reduce to about 5% once the lead- 



ing logarithmic ISR corrections are included in RacoonWW (RacoonWW IBA). The 
remaining differences should be due to the still quite different treatment of visible photon 
radiation in RacoonWW and YFSWW3: in contrast to RacoonWW, YFSWW does not 
include the complete lowest-order matrix elements for e"'"e~ ^fl- Instead, the photon 



radiation from the final state is treated via PHOTOS [|28[. In particular, for small photon 
energies, where the differences are largest, the non-factorizable contributions, which are 
not yet included in YFSWW3, might play a role. 

In Figures |8|-p!0| we extend the comparison of the photonic distributions between 
YFSWW3 and RacoonWW to 500 GeV. Here the difference is typically at the level of 
10% and in general not reduced by the inclusion of ISR for e"^e~ — > 4/7 in RacoonWW, 
i.e. the agreement without ISR in RacoonWW was accidentally good. One should also 
recall that the diagrams without two resonant W bosons (background diagrams) become 
more and more important at higher energies. Thus, the increasing difference between 
YFSWW3 and RacoonWW for higher energies could be due to a less efficient descrip- 
tion of final-state radiation by the effective treatment with PHOTOS. 
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Figure 6: Distribution in the cosine of the polar angle of the photon w.r.t. the e"*" beam 
for the process e"'"e~ ud/x~z/^7 at a/s = 200 GeV 
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Figure 7: Distribution in the angle between the photon and the nearest charged final-state 
fermion for the process e"'"e~ ud/i~z/^7 at y/s = 200 GeV 
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Figure 9: Distribution in the cosine of the polar angle of the photon w.r.t. the e"*" beam 
for the process e"'"e~ — > ud/i~z/^7 at y/s = 500 GeV 
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Figure 10: Distribution in the angle between the photon and the nearest charged final- 
state fermion for the process e+e" — *• ud/x~i>^7 at y/s = 500 GeV 



5.2 Standard Model predictions 

We now discuss the predictions of RacoonWW for various observables in the SM. 
Here and in the following, the width is always calculated from the input parameters in 
lowest order in the scheme including naive QCD corrections (Fw = 2.09436 GeV). 
Naive QCD corrections are included in all results, in particular, also in the Born results. 

The results in the LEP2 energy range were obtained with the ADLO/TH cuts as 
defined in Ref . [ || , those at y/s = 500 GeV with the cuts 

^(Z, beam) > 10°, 9{l,l') > 5°, 9{l,q) > 5°, 

^(7,beam) > 1°, ^(7,/)>5°, ^(7,g)>5°, 

> 0.1 GeV, El > 1 GeV, Eg > 3 GeV, 

m(g,g') > O.lGeV, ^(g, beam) > 5°, (5.1) 

where 9{i,j) specifies the angle between the particles i and j in the LAB system, and 
/, q, 7, and "beam" denote charged final-state leptons, quarks, photons, and the beam 
electrons or positrons, respectively. The invariant mass of a quark pair qq' is denoted by 
m{q, q'). The cuts ( |5.1| ) differ from the ADLO/TH cuts only in the looser cut on m(g, q') 
and in the additional cut on ^(g, beam). 

In Figure ITT] we present the total cross section for the process e"'"e~ ud/i~z/^7 in 



the LEP2 energy range. On the left-hand side we show the absolute prediction in lowest 
order (Born), including ISR (ISR), and including in addition the Coulomb singularity 
according to variant 1) (IB A) discussed in Section ^ On the right-hand side we give the 
corrections relative to the lowest order including in addition a curve with the Coulomb 
singularity according to variant 2) (IBA2). The Coulomb singularity reaches about 5% 
at threshold and decreases with increasing energy. The effect is comparable to the one 
for the process without photon. The two variants for the implementation of the Coulomb 



14 




150 160 170 180 190 200 210 
V^[GeV] 




150 160 170 180 190 200 210 
V^[GeV] 



Figure 11: Total cross section for the process e^e ud/i as a function of the CM 
energy. 



singularity show hardly any difference. Consequently, we will always use variant 1) in the 
following. 

In Figures p!2|-|l^ we present the distributions in the photon energy E^, in the cosine 
of the polar angle 6^ of the photon w.r.t. the e^ beam, and in the angle 6^f between 
the photon and the nearest charged final-state fermion for ^/s = 200 GeV. The left-hand 
sides contain the absolute prediction for the process e"'"e~ udyU~z/^7 in lowest order 
(Born) and including the ISR corrections and the Coulomb singularity (IBA), and for 
the process e"'"e~ — > ude~z/e7 including these corrections. The relative corrections (right- 
hand sides) are typically of the order of —10% wherever the cross sections are sizeable. 
Relative to the corresponding lowest-order results, the corrections to e"'"e~ ude~i>e7 
would practically be indistinguishable from the relative corrections to e"'"e~ ud/i~z/^7. 
We therefore prefer to plot the corrections to e"'"e~ — ude~t'e7 normalized to the lowest- 
order of e"'"e~ ud/i^z/^7 in order to visualize the effect of the "background" diagrams 
contained in e^e^ —>■ ude~t'e7- As can be seen, this effect is comparable to the radiative 
corrections but of opposite sign. 

In Figures p!5|-[l7| we show results for ^/s = 500 GeV. Here, the left-hand sides contain 
the absolute prediction for the processes e+e" — >• ud/i~z/^7 and e"'"e~ ude~t'e7 in low- 
est order (Born) and including the ISR corrections and the Coulomb singularity (IBA). 
Note that here the distributions differ sizeably between the two processes. Therefore, on 
the right-hand sides, the IBA predictions for both processes are normalized to the corre- 
sponding lowest-order predictions. Where the cross sections are sizeable, the corrections 
are about +10% for e+e~ ud/i~z/^7 and +5% for e"^e" ude~i/e7- They are larger 
where the cross sections are small. 
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Figure 14: Distribution in the angle between the photon and the nearest charged final- 
state fermion for the processes e"'"e~ — > ud/i~z/^7 and e"'"e~ — > ude'i^eT at ^/s = 200 GeV 
("Born(/x)" indicates that the Born cross section for e"'"e~ —>■ ud/i~z/^7 is taken.) 
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Figure 16: Distribution in the cosine of the polar angle of the photon w.r.t. the e"*" beam 
for the processes e+e~ — > ud/^~t'^7 and e"'"e~ ude~z/e7 at i/i = 500 GeV 




20 40 60 80 100 120 20 40 60 80 100 120 



e^f [deg] e^f [deg] 

Figure 17: Distribution in the angle between the photon and the nearest charged final- 
state fermion for the processes e"'"e~ — ^ udyU~z/^7 and e"'"e~ — >• ude~z/e7 at y/s = 500 GeV 
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5.3 Predictions with anomalous quartic couplings 



Since the matrix element depends hnearly on the anomalous quartic couplings a^, the 
cross section is a quadratic form in the Oj. Therefore, it is sufficient to evaluate the cross 
section for a finite set of sample values of the anomalous quartic couplings in order to 
get the cross section for arbitrary values of these couplings. We restrict ourselves here to 
the semileptonic process e+e~ — > ud/i~i/^7 and include ISR and the Coulomb singularity 
(variant 1). We use the cuts 

E^>5GeV, I cos I < 0.95, | cos ^^/| < 0.90, m(/, /') = Mw ± 2rw, 

(5.2) 

where is the energy of the photon, 6^ the angle between the photon and the beam axis, 
9^f the angle between the photon and any charged final-state fermion /, and m(/, /') the 
invariant mass of the fermion-antifermion pairs that result from W decay. In the compu- 
tation of m(/i, Vy) the momentum of the neutrino is set equal to the missing momentum, 
since the neutrino is not detected, i.e. the energy loss in the ISR convolution ( [4.2|) is 
implicitly included in the neutrino momentum. 

We first study the infiuence of the AQGC oq, Cc, Cn, ao, and on the cross section 
at \/s = 200 GeV and 500 GeV separately. Figure |1^ shows the cross section normalized 
to the SM value as a function of each of these couplings for all the other Cj's equal to 
zero. The asymmetry results from the interference between the SM matrix element and 
the matrix element of the AQGC, which is suppressed for the CP-violating couplings On 
and oq. The asymmetry is small for oq and fin and only visible at ^/s = 500 GeV for in 
Figure |18|, but sizeable for Oc.^ The cross section is most sensitive to and fio and least 



sensitive to and fin. 

In order to illustrate the potential of LEP2 and a linear e"'"e~ collider in putting limits 
on the AQGC, we consider the following two scenarios: an integrated luminosity C = 
320 pb"^ at y/s = 200 GeV and C = 50 fb"^ at = 500 GeV. The corresponding total 
SM cross sections to e"'"e~ uciyU^z?^7 are 16.69 fb and 7.64 fb, respectively. Assuming 
that the measured number of events is given by the SM cross section asM = <^{cLi = 0) 
and the experimental errors by the corresponding square-root, we define 



2_ iN{a,)-Ny (a{a,) 




X' = = - 1 (TsuC, (5.3) 

where N{ai) is the number of events that result from the cross section with anomalous 
couplings. Since the square- root of this distribution is a quadratic form in the Oj, 
the hypersurfaces of constant form ellipsoids. The la limits resulting from = 1 
on individual couplings can be illustrated by projecting the ellipsoids into the planes 
corresponding to pairs of couplings. Instead of the projections, often the sections of the 
planes with the ellipsoids are used. Note that the ellipses resulting from projections are 
in general larger and include those ellipses resulting from sections of the planes with the 
ellipsoids. Since the correlations are small for the cases under consideration, the difference 



"^The sign of the asymmetry differs from the results of Ref. [|^, since the couphngs oq and Oc have 
been implemented [ ^ in EEWWG with a sign opposite to the definitions in Refs. [|7[ 0, which agree 
with our choice. 
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Figure 18: Impact of the AQGC oq, etc, Qn, Sq, and an on the cross section for e"''e^ 
ud/i~z/^7 at a/s = 200 GeV and 500 GeV. Only one of the AQGC is varied while the 
others are kept zero. 
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between both types of ellipses is also small. In the following figures we include both the 
projections and the sections of the ellipsoids using the same type of lines. 

In Figures [1^ and ^ we show some la contours for various pairs of Oj. In addition 
we list the la limits derived from projecting the ellipsoids. Since the effects of and do 
and of On and d^ on the cross section are equal up to relatively small interference terms, 
also the corresponding contours of these couplings with other couplings are of similar 
size. For transparency we omitted some contours involving a^; for 200 GeV (Figure |T^) 
these contours practically coincide with the ones for d^, for 500 GeV (Figure 12^) the 
contours for On are of the same size and shape as the ones for d-^ but shifted to become 
approximately symmetric w.r.t. —dn- The best limits can be obtained for oq and 

oq. The correlations between the different couplings are in general small, and only oq 
and ttc show a noticeable correlation. The limits obtainable at a linear collider are by 
about a factor of 200 better than those obtainable at LEP2. This improvement reflects 
the enhanced sensitivity of the cross section on the anomalous couplings at high energies, 
which can also be seen in Figure [l^, and to a smaller part the higher luminosity. 



6 Summary 

We have calculated all lowest-order amplitudes for e+e~ —>■ 4/7 with five different gen- 
uine anomalous quartic gauge-boson couplings that are allowed by electromagnetic gauge 
invariance and the custodial SU(2)c symmetry. These couplings include the three opera- 
tors Co, Cc, and £n, which have been constrained by the LEP collaborations by analysing 
WW7 production, and two additional P-violating couplings, one of which conserves CP. 
The five anomalous couplings have been incorporated in the 4/(7) Monte Carlo generator 
RacoonWW. We have calculated the dependence of the cross section for e"''e~ 4/7 
on the anomalous quartic couplings and illustrated the typical size of the hmits that can 
be obtained for these couplings at LEP2 and a 500 GeV e"'"e~ collider. 

Moreover, we have implemented the dominant leading electroweak corrections to 
e"'"e~ — ^ 4/7 into RacoonWW. These include initial-state radiation, the dominant uni- 
versal effects originating from the running of the couplings, and the Coulomb singularity 
for processes involving W-boson pairs. We have compared the corresponding predictions 
with existing calculations, as far as possible, and investigated the numerical impact of the 
dominant corrections. 

With the additions described in this paper, RacoonWW is a state-of-the-art Monte 
Carlo generator for the classes of e"'"e~ 4/ and e"'"e~ 4/7 processes with arbitrary 
massless four-fermion final states, both for the Standard Model and including anomalous 
quartic gauge-boson couplings. 
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Appendix 

A Some corrections to the generic construction of e+e" 



4/7 amplitudes 



In Ref. [ 1^ we have constructed the amplitudes for all e"^e~ 4/7 reactions from the 
two basic channels CCa and NCa, which are also specified in Section ^ Here we take the 
opportunity to correct two mistakes in the corresponding formulas: 

• Equation (2.24) of Ref. [ ^ is only correct for down-type fermions /, while some 
arguments have to be interchanged for up-type fermions. The correct formula is 

-^Cc7NCb''''''''''''(P+'P-' ^1' ^2, h) 



^ k3,-ki,ki,-p^, 
{p+,p^,ki,k2,k3,k4) 



for Ilj 



-1/2, 



» y Cr+ ,0-_ ,(71 ,(T2 ,0-3 ,0-4 



for II f = +1/2. 



(A.l) 



The error affected the evaluation of the final states i^ei^ei^ni^fi and z/gt'eUU in Table 1 
of Ref. [|^ at the level of 0.2-0.4%. The corrected results for Table 1 are 



cr/fb 


e+e- 4/ 
running width 


e+e- ^ 4/ 
constant width 


e"'"e 4/7 
constant width 




8.339(2) 


8.321(2) 


1.511(1) 


t'e^'eU U 


23.91(2) 


23.90(2) 


6.79(3) 



For the final states i^e^e^fi^^^il and UeiJeUn'j no change is visible in the numerical 
results within the integration errors after the correction. The numerical smallness 
of the correction is due to the fact that the two cases in ( A.l ) differ only in the 
contribution of a non-resonant background diagram which is suppressed. 

Equation (2.25) of Ref. [ Q contains some misprints. The correct formula is 
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However, the numerical evaluations for the corresponding z/eZ/eZ/eZ/e(7) final states 
were based on this correct form. 
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